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Abstract. We undertake an exact numerical study of the effects of disorder on the Anderson localization 
of electronic states in graphene. Analyzing the scaling behaviors of inverse participation ratio and geo- 
metrically averaged density of states, we find that Anderson metal-insulator transition can be introduced 
by the presence of quenched random disorder. In contrast with the conventional picture of localization, 
four mobility edges can be observed for the honeycomb lattice with specific disorder strength and impurity 
concentration. Considering the screening effects of interactions on disorder potentials, the experimental 
findings of the scale enlarges of puddles can be explained by reviewing the effects of both interactions and 
disorder. 

1 Introduction with weak long-range impurities, but the intervalley scat 

taring and chirality breaking scattering should be consid 
As a two-dimensional (2D) allotrope of carbon [F5], grapheneered in the presence of strong short-range disorder. There 



' holds great pronrise for replacing conventional senricon- 
. ductors on account of its unique electronic properties. In 

■ graphene, the cr— band formed by the sp^ hybridized or- 
bitals determines the robustness of the honeyconrb struc- 

' ture, and the half-filled 7r-band provided with the Dirac- 
' like electronic excitations is responsible for the unusual 
, transport properties [3|. It has been observed in some ex- 

■ periments [4 5 6 7 8 9J that disorder has remarkable effects 
[ on the unusual electronic properties of graphene, and thus 
' there has been a great deal of interest in recent years con- 
cerning the role of disorder in graphene. 

Some experiments have confirmed the existence of metal- 
' insulator transition (MIT) in disordered graphene |4I5I6) . 
, which suggests that the one-parameter scaling theory [TU] 

■ may break down [llj. In accordance with the prediction 
I of Klein paradox [12J, the Dirac fermions are found to 
. be robust against disorder in the 2D single valley Dirac 

nrodel. In addition, the one-parameter scaling in single- 
valley Dirac Hamiltonian shows that |13I14| , the logarith- 
' mic derivative dlna/dhiL — (3(cr) is a positive function 
' and the scaling flow has no fixed point, indicating that 
the system always scales towards a metal. This scenario 
is different from what would be expected for the conven- 
tional 2D electron systems, where all states are localized 
for arbitrary weak disorder |10| . 

Recently, some nunrerical methods have been adopted 
in the study of the massless Dirac model with random 
potential [14 15 16 17 18J or the Anderson tight-binding 
model |19I20I21I2 2 23 24 25I26I27J . It is already clear that 
the single valley Dirac approximation is valid for graphene 
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fore, the full tight-binding description of graphene with 
Anderson disorder becomes a popular alternative model. 
Some very recent simulations |22I23I24] have provided new 
evidence for the existence of MIT in graphene with short- 
range disorder by obtaining the mobility edges near the 
Fermi energy. However, a contrary view has been sup- 
ported by the transfer-matrix approach |25| and kernel 
polynomial method [27] that all states in the Anderson 
tight-binding nrodel are localized for arbitrary weak dis- 
order. The possible reason for the disagreement to arise 
is that different measures and scaling rules are adopted 
in the above studies in an effort to distinguish localized 
states from delocalized ones. 

In this paper, we introduce a new approach to scale 
the inverse participation ratio (IPR), which overcomes the 
unstable and unreliable problems caused by the primitive 
scaling of IPR with the negative first or secoird power of 
lattice size. We study the tight-binding Hamiltonian of the 
honeycomb lattice with randomly fluctuated disorder po- 
tentials, and an unique picture of Anderson MIT has been 
obtained. As expected, we flnd that all electronic states in 
graphene are localized by strong disorder. However, four 
mobility edges can be observed when we decrease the dis- 
order strength until under a critical value Wc, and the elec- 
tronic states within two nearest neighbor mobility edges 
are entirely extended or localized alternatively. Therefore, 
MIT can be introduced by changing the carrier density 
to move the Fermi surface across a certain mobility edge. 
Our results are in agreenrent with the very recent exper- 
imental suggestions, where a transition from localization 
to antilocalization has been achieved by decreasing the 
carrier density in graphene sanrples |5]. 
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On the other hand, the sources of the coexistence of 
electron and hole puddles observed by the scannable single- 
electron transistor (SET) have been proved to be both the 
substrate-induced structural distortion as well as chemi- 
cal doping. While the length scale of density fluctuations, 
which is more than 150nm [7], is extremely larger than 
the disorder length scale introduced by the short-range 
scatters. In this paper, we show that the unexpected large 
disorder length scale observed by SET can be explained 
by the screening effect of interactions on disorder, suggest- 
ing that the electron-electron interactions should be taken 
into account. Besides, the delocalization effect of interac- 
tions through screening disorder potentials can also intro- 
duce the zero-bias anomaly at Fermi surface [28129130] . 

In this paper we address the following issues related to 
the effect of disorder in graphene: (1) Whether there exists 
the mobility edge to separate the localized states from the 
delocalized states in disordered graphene? A distinctive 
scenario of Anderson MIT with four mobility edges have 
been observed in graphene with quenched random disor- 
der; (2) How to measure the localization length of a local- 
ized electronic state by IPR? We present a new method to 
set the lattice size scaling of IPR, which can give the local- 
ization length of electronic state precisely; (3) What is the 
role of interactions in disordered graphene? The screening 
effect of the short-rang interactions on disorder potentials 
is found to have strong affection for the scale enlarges of 
puddles. 

The paper is organized as follows: In section. 12.11 we 
present the fully tight-binding description of the disor- 
dered graphene. In section. 12.21 and 12.31 we discuss the 
lattice size scaling of the IPR and the geometrically aver- 
aged density of states (GADOS) respectively, and show 
the scaling method relevant for the subsequent discus- 
sions. Next we show our findings regarding the static dis- 
order in graphene in secrion. |31 In addition, the screening 
effects of the electron-electron interactions on the disor- 
der potentials are presented in section. SI The principal 
findings of this paper are summarized in section. [S] 



2 Theory 

2.1 The Anderson Tight-banding Model 



The Anderson tight-banding Hamiltonian [31132] for the 
disordered graphene can be expressed as 



Ham = -< 51 ^' "^J + Xl ^«^»^'^*' 



(1) 



where Ci {c[) are annihilation (creation) operators of a 
fermion at site i, ti represent the on-site disorder ener- 
gies, and t denote the hoping integrals between nearest 
neighbor (NN) sites, which are set to be unit in our calcu- 
lations. We take the exact numerical approach to find all 
eigenenergies and eigenfunctions of the Hamiltonian given 
by Eq. ([T|) for a finite hexagonal lattice with the periodic 
boundary conditions. Besides, the numerical calculations 



are performed with regard to the following three different 
disordered cases. 

Firstly, in the presence of quenched random disorder, 
we introduce the box distributed Anderson disorder, and 
the on-site disorder energies e^ in Eq. ^ are assumed to be 
independent random variables distributed between —W/2 
and W/2. Therefore, the probability distribution of e^ is 
given by 



P(e) = W~^0{W/2-\t\), 



(2) 



and the disorder strength is parameterized by the width 
W. 

Secondly, to take into account the infiuences of the im- 
purity potentials to the adjacent sites [5T], the correlation 
length of disorder C is introduced and the on-site disorder 
energies ti in Eq. ([T]) can be replaced by the correlated 
disorder energy e^ as 



N 



e, =^Qea;p(-|r, -r,|VC^). 



(3) 



1=1 



Where e; represent the "uncorrelated" disorder energies of 
the impurities at sites I, which have been defined in the 
first case. 

Finally, in order to generalize this Hamiltonian by in- 
cluding electron-electron correlations, one must add an 
interaction term Hint onto the Anderson tight-binding 
model: 

H = HAM+H„,t, (4) 



with 



Hj, 



-^E 



ni^riii 



E^^ 



ijTl-lTlj . 



(5) 



Here ni^r — cj^Ci^ is the local charge-density operator with 
a =t, J,, U represent the on-site interactions, and Vij de- 
note the non-local interactions between electrons at sites 
i and j. To compare with the experimental finding about 
the enlarges of puddles [7], here we consider a case of the 
binary alloy disorder, in which the probability distribution 
function of the on-site disorder energies is given by 



p(e,) = xS{e, -W) + {1- x)S{e.,), 



(6) 



where x is the fraction of the lattice sites with energies 
Ei — W, and W represents the on-site energy splitting. 



2.2 The Lattice Size Scaling of IPR 

One of the widely used quantities to measure the Anderson 
localization of electronic states in disordered systems is 
called IPR, which is defined as [33I34| 



N 

El 



/2(i^„) = V|tf'„(r.)l^ 



(7) 



where En (n=l,...,N) and ^nifi) are the eigenenergies and 
eigenfunctions of a disordered finite lattice with N sites. 
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Since the eigenenergies of a finite lattice are always dis- 
crete, it is more convenience to introduce a continuous 
energy-dependent IPR as 



^N 



l2iL0,N) 



EtiOi; 



-En\)l2{En) 



EL0(i 



w- 



En\) 



(8) 



where uj represents energy, and 7 is a very small energy 
bin for the average. Since IPR of extended states are pro- 
portion to 1/ L'^ and go to zero in the infinite lattice limit 
L -^ 00, it is not difficult to identify extended states by 
doing lattice size scaling of IPR. Here d denotes the di- 
mension of the system, and L represents the lattice size 
with N — L"^. However, the lattice size scaling of IPR for 
a localized state is found to be more complicated than 
expected |3KI;^6l37l:i8| . 

The wave function of Anderson localized states de- 
cays exponentially as |V'('")| ^ exp{—\r — ro|/^;oc), where 
^joc is the localization length. When the lattice size L is 
much larger than the localization length, IPR is a size- 
independent constant J", and the localization length can 
be obtained easily by S^ioc — (^2)^^^'* L^SJ- However, we 
can not always perform calculations satisfying the condi- 
tion L ^ ^loc on account of the limited capacity of the 
numerical calculations. When the lattice sizes are compa- 
rable to or even smaller than the localization length S^ioc, L 
becomes a variable to the function of IPR. Later a simple 
lattice size scaling formula 1/L" has been suggested. Un- 
fortunately, this IPR scaling can mistakenly regard a local- 
ized state as an extended state in some cases. For example, 
the exact scaling formula of IPR in one-dimensional (ID) 
systems with box distributed disorder has been obtained 
as ESI 



12(0;, L) =/2(w,oo)coth(L/2^;oc(w)), 



(9) 



with S,ioc{^) = I2 (w,c»). It is obvious that the simple 
formula 1/L" is not a proper approximation in consider- 
able large lattice size region to scale IPR. Certainly, there 
are also exceptions for the conditions with L 3> ^loc and 
L <C ^loci where we can choice a = or a = 2 respectively. 
Unlike the ID disordered system, the exact expression 
for IPR has not been obtained for the 2D finite systems. 
Therefore, we introduce a Taylor series to scale IPR by 



12(0;, L) = /2(w,oo)-|- 



ai{Lo) 02(0;) 03(0;) 



L 



L2 



L3 



(10) 



where a„(a;) is the n-th Taylor parameter. The minimum 
radius of convergence is found to be Rmin = 0.1 in our 
studies, and to make all Taylor series adopted conver- 
gent, we have to do calculations when lattices meet the 
condition of L > 1/ Rmin — 10. For convenience, we em- 
ploy a polynomial formula consist of the front five terms 
in Eq. p^ since the contributions of high-order terms 
are negligible. As shown in Fig. [Tl the new fitting is rea- 
sonable and pratical since the intercept of an extended 
state is found to be zero or even a very small negative 
number, whereas the fitting curve of a localized state has 
a positive intercept in the infinite-lattice limit. We have 



3" 




l/L 

Fig. 1. (color online) The lattice size scaling of IPR 
l2{uj,L) for a localized (circles) and an extended (trian- 
gles) state. The black lines are the fitting curves obtained 
by the polynomial formula shown in Eq. (|10p . and the 
red lines show the extensions of the fitting curves to the 
infinite lattice limit L -^ 00. 



checked the above fitting method by the disordered cubic 
lattice [39], and good agreements have been achieved with 
the accepted scaling theories [T^. Especially, we find that 
the fitting curves of extended states can be approximated 
by l2{uj,L) oc 1/L^, in according with the prediction of 
some other theories. On the other hand, when the elec- 
tronic state is localized, the five parameters are all found 
to have non-zero values, but the low-order terms play the 
main roles. As show by the solid line in Fig. [U there ex- 
ists a finite intercept 12(0;, 00) for a localized state, and 
its localization length can be obtained by 



^loci(^) 



^7/2(^,00) 



(11) 



To sum up, the new scaling method of IPR has the 
advantage in distinguishing explicitly the localized states 
with extended states, and in addition we can acquire the 
localization length precisely by the intercept obtained in 
the infinite-size limit. 



2.3 Measure Anderson Localization by GADOS 

It has been proposed that Anderson localization can be 
measured by an order parameter Pg{uj), called GADOS, 
which is obtained by geometrically averaging the local 
DOS (LDOS) of each site as 
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(13) 



Where p{ri,E„) = \>Fn{ri)\'^ is the LDOS at site i for the 
n-th eigenstate tf'n, Ns represents the number of disorder 
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Fig. 2. (color online) The energy dependencies of the lo- 
calization length S,ioc(j-^) of graphene with Anderson dis- 
order. The box distribution is adopted for the on-site en- 
ergies of every sites with disorder strength 14^=2.0 (red 
triangles), 3.0 (blue squares), and 4.0 (black circles). The 
dotted lines denote the mobility edges. 
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Fig. 3. (color online) The lattice size scaling of GADOS 
of electronic states with different energies: (a) uj = 3.175 
(black circles), and (b) oj =0.5 (black triangles) and 2.5 
(blue squares). Insert: the scaling of GADOS close to the 
Dirac point with w = 0.0125 (green circles). The disorder 
strength of Anderson disorder is W=2.0. 



configurations to be averaged, and the energy bin 7 plays 
the same role as in Sec. 12.21 

For the infinite-dimensional disordered system, it has 
been defined that Anderson transition happens when pg (ep) 
vanishes completely at Fermi surface ep- This criterion 
has been introduced to investigate the competition be- 
tween Anderson transition and Mott MIT in the infinite- 
dimensional systems within the dynamical mean-field the- 
ory |40I41| . While, in the numerical calculation for a finite 
lattice of low-dimensional disordered system, a finite en- 
ergy bin 7 has to be introduced. As a result, GADOS is a 
function of lattice size, and its lattice size scaling should 
be introduced to detect the Anderson MIT. 

In Ref. [3S], we have examined the lattice size scale of 
GADOS and found that, for any nonzero energy bin 7, 
Pg decays exponentially with the increasing of lattice size 
for a localized electronic state. On the contrary, there is 
no significant variation of pg with the increasing of lattice 
size for a delocalized state. Therefore, it is not difficult to 
distinguish localized states from extended states by the 
lattice size scaling of GADOS. In this paper, we use GA- 
DOS scaling to check our results obtained by IPR scaling 
about whether the electronic states are localized or not, 
and good agreements have been achieved. 



3 Localization in disordered graphene 

The scaling theories of localization have proved that all 
electronic states should be localized in the conventional 
2D systems with arbitrary weak disorder ^U\. Compar- 
ing with the parabolic DOS of conventional 2D system, 
graphene has a special band structure with zero weight 
at Dirac points. Applying the Dirac model with random 
potentials to the weakly disordered graphene, some stud- 
ies have confirmed the existence of the antilocalizatized 
electronic states [15 14J because of the insensitivity of the 
Dirac fermions to the disordered external electrostatic po- 
tentials. On the other hand, strong disorder can mani- 
festly enhance the DOS at Dirac points and strongly affect 



the electronic properties. For this reason, our study will 
concentrate firstly on the evolutions of electronic states 
in the whole energy band with the increasing of disorder 
strength. 

We study the Anderson tight-binding model on a hexag- 
onal lattice by exact numerical simulations, where the 
largest lattices could be 120x120. Applying the IPR scal- 
ing method introduced in Sec. l2.2[ the localization lengths 
of localized electronic states are obtained for different dis- 
order strength W. As shown in Fig. [51 it is obvious that 
disordered graphene has a quite different scenario of An- 
derson localization than the conventional 2D systems with 
disorder. When disordered strength ly = 2 or 3, we find 
energy regions consist of the localized electronic states 
separated by the mobility edges from the regions of ex- 
tended ones, where the physical observable is the finite lo- 
calization lengths of localized states. Close to the top and 
bottom of the whole energy band, we obtain two symmet- 
rical mobility edges at a; = ±Wci • In addition, two mo- 
bility edges appear near the Fermi surface {ep = 0) with 
UJ = ±a;c2(wc2 < Wci), suggesting the electronic states 
around the Fermi surface are also easier to be localized by 
the scattering with the disorder potentials. Furthermore, 
all electronic states can be localized when the disorder 
strength is increased larger than W^ = 3.10. 

The region close to the Dirac points is very important, 
and people show great interests in the possible existence 
of the antilocalized electronic states at there. However, 
the scaling results are affected magnificently by the rela- 
tive big errors for IPR since the DOS is quite low in the 
vicinity of Dirac points. In addition, it is well known that 
to distinguish an extended state from a weakly localized 
state with large localization length is very difficult. To 
make our finding more convincing, we also use the scaling 
of GADOS to measure Anderson localization in disordered 
graphene. We find very good agreements between the scal- 
ings of GADOS and IPR. 

In Fig. [31 we show the lattice size dependencies of GA- 
DOS for some representative energies when W = 2. In 
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Fig. 4. (color online) (a) The energy dependencies of the 
critical disorder strength Wc{uj) of graphene with corre- 
lated disorder; (b) The effects of the correlations of disor- 
der on DOS around Fermi surface w = 0. The solid (black), 
dashed (blue) , and dashed-dotted (red) lines represent re- 
spectively the results of point defects, NN scatters and 
both NN and NNN scatters. 



Fig. [3l^a), GADOS of w = 3.25 increases manifestly with 
the increasing of 1/L, suggesting the state is localized. On 
the contrary, GADOS of the other two extended states 
with energy w = 0.5 and 2.5 are both lattice size indepen- 
dent as shown in Fig.[3l^b). When w = 0.0125, we have per- 
formed the scaling of GADOS very carefully and precisely 
by averaging more than two hundreds of disorder config- 
urations for lattice with L = 120. As shown by the insert 
in Fig.[2I^b), the lattice size independent behavior of GA- 
DOS indicates clearly that the electronic states quite close 
to the Dirac points are still extended, in accordance with 
the result obtained by the IPR scaling. Our result sug- 
gests that the prediction of single valley model, i.e. Dirac 
Fermions can not be localized by disorder, is reasonable in 
weakly disordered graphene since no intervalley scattering 
can be stimulated by very weak disorder. In addition, it is 
obvious that the Anderson MIT can be introduced by only 
changing the carrier density to move the Fermi surface 
from extended region to localized region with fixed disor- 
der strength. Recently, a transition between localization 
and antilocalization has been observed experimentally in 
graphene when the carrier density is decreased [B] . There- 
fore, our result about Anderson localization with four mo- 
bility edges help clarify this experimental observation of 
MIT in graphene. In addition, our findings suggest that 
the delocalization is closely related to the linear dispersion 
near the Dirac points which may suppress the inter- valley 
scattering. And further investigations are to be done in 
the future study. 

The electrostatic coulomb potentials of the surface ab- 
sorptions or the adatoms in the substrate for graphene 
should be described by Eq. ([3]), in which the correlations 
between the disorder potentials of individual sites are con- 
sidered. Owing to the screening effect of electrons, the 
long-range correlations of disorder could be neglected. Here 
we name the uncorrelated and short-range correlated dis- 
order as point defects and short-range scatters respec- 
tively. In Fig.m we show the effects of the nearest-neighbor 
(NN) and next NN (NNN) scatters on the locaHzation of 
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Fig. 5. (color online) The evolutions of DOS of graphene 
with disorder strength W = A and short-range interac- 
tions: (a) there are only on-site interactions U =0.5 (black 
dashed line), 1.0 (red dashed line), and 2.0 (blue solid 
line), and (b) the on-site interactions are fixed asU — 2.0, 
and the NN interactions are given as Vi = 0.0 (black 
dashed line), 0.5 (red dashed line), and 1.0 (blue solid 
line) . 



electronic states. The correlation length C is set as 1.5a, 
which is the same as the parameters chosen in Ref. |42| . 
and the blue dashed and red dash-dotted lines reveal the 
energy dependence of the critical disorder strength Wduj) 
for graphene with the NN scatters only and both NN and 
NNN scatters, respectively. For the convenience of com- 
parison, we also show the results of uncorrelated point 
defects (black solid line). It is obvious that the picture of 
Anderson MIT for short-range scatters is the same as that 
of the point defects. While, the short-range correlations 
of disorder have strong effects to localize the electronic 
states. As to why the critical strengthes of disorder Wc for 
localizing the states are greater for uncorrelated disorder 
than the correlated ones, the effective disorder strength is 
actually enhanced a lot as we add the influences of the 
impurity potentials of the adjacent sites onto the uncorre- 
lated on-site energy of a certain site directly. As a result, 
the localization of electronic states is enhanced accord- 
ingly when the re-scaling is not introduced to reduce the 
original disorder strength. 



4 The interactions and electron-hole puddles 
in disordered graphene 

In pure graphene, the influence of short-range interactions 
on electronic properties has been investigated by the dy- 
namical mean-field theory (DMFT) [33]. The Dirac sea 
state is found to remain stable against local many-body 
interactions since the interactions in graphene are much 
smaller than the Mott critical value Uc = 13. 3t [33|- In the 
present of disorder, the combine effects of interactions and 
disorder will have strong influences on the LDOS and also 
the localization of electronic states. It has been proved 
recently by the statistic DMFT calculations [3S] that the 
Hartree-Fock approximation (HFA) could give reasonable 
results for the conventional 2D Anderson-Hubbard model 
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when the interactions U are smaller than the energy band- 
width D. Based on the predictions for the many-body ef- 
fect in graphene given by Ref. [3_, we study the Anderson- 
Hubbard mode with weak interactions in the region of 
U = 0.5t to AM within HFA. 

In conventional 2D disordered systems, the zero-bias 
anomaly at Fermi surface ^28.29 30J arises from the inter- 
play between disorder and interactions, which indicates 
the delocalization effects of interactions on localized elec- 
tronic states 131]. In the same manner, we observe the 
DOS of disordered graphene at Fermi surface to show 
the delocalization effects of the short-range interactions. 
Firstly, we only consider the effects of the on-site interac- 
tions U, and as shown in Fig. [SJa), significant decreasing 
of the DOS around the Fermi surface appears when we 
increase U from 0.5 to 4.0. On the contrary, the long-rang 
interactions are found to have less effect to the zero-bias 
anomaly at Fermi surface. In Fig. [5jb), we plot the DOS 
for graphene with different NN interactions and fixed U, 
and it is obvious that even the effects of NN interactions 
are still weak and negligible. Therefore, we predict that 
the on-site interaction is the key point in the study of the 
screening effects of interactions in disordered graphene. 
The local effect of interactions can be understood by calcu- 
lating the screened potential Vi = ei{l — Uxu), where the 
local charge susceptibility is defined as xa — —dni/det. 
It is obvious that the screening effect in localized phase 
is expected to be less than in metallic phase since xu is 
restrained for site with small LDOS at Fermi level. We 
find in the paramagnetic phase that the on-site interac- 
tions have strong effects to delocalize the electronic states 
by increasing the localization lengths accordingly. 

Since The LDOS at Fermi surface ep is detectable, 
to compare with the experimental results, we calculate 
the evolution of piri^ep) with the on-site interactions U 
for a particular binary disordered configuration. As shown 
in Fig. [6] and Fig. [7l the length scale of density varia- 
tions of LDOS for U — 0.0 is much smaller than that 
for U = 2.0 and 4.0. Compared with the quarter-filling 
case, the delocalization effects of interactions are not re- 
markable when the system is half-filled. It is obvious that 
the states at Dirac points are strongly localized at half- 
filling with W = 4:. This further proves the results we 
obtained in the previous section from a different perspec- 
tive. The regions of electron-rich and hole-rich puddles 
have been observed by the experiment of using scanning 
single-electron transistor to map the LDOS of graphene 
sheet ^ . The smallest length scale on which density vari- 
ations is observed roughly 150nm, which is apparently sig- 
nificantly larger than the intrinsic disorder length scale 
as approximately 30nm in the graphene samples. There 
arises the question of why there exists such a big differ- 
ence between the length scales of the LDOS and intrinsic 
impurities. Our results predict that the interactions play 
an important role in the electron-hole puddles observed 
by experiment, where the screening effects of the on-site 
interactions on disorder potentials can enlarge the length 
scale of LDOS significantly. 



00 - 


^^^1 


50 - 


^H^ 


- 


^^H 



100 



50 



(a) 



'■ > M 


w 



50 

X 



100 ^ 0.0 

^~ 4.0 



(b) 



50 

X 



100 



I 0.0 
. .002 
I .004 
1 .006 

I .008 



00 - 
50 - 


n^^3 


1 


1 


li 






1 




a 








50 


100 







(c) 




X 





00 - 


'W^ 




50 - 


. » i*. A.' ^ ' 1 


f ' 


- 





0.0 

.002 
.004 
.006 



(d) 



50 

X 




Fig. 6. (color online) (a) The disorder configuration of 
an 108x108 lattice with binary disorder. The on-site en- 
ergy splitting is ly = 4; The corresponding LDOS of an 
eigenstate right at the Fermi level for different on-site in- 
teractions: (b) U = 0.0, (c) 2, and (d) 4.0. The systems 
are all half-filled with impurity concentration x — 0.2. 



5 Conclusions 

We have studied numerically the Anderson tight-banding 
model of the finite hexagonal lattices and found that, un- 
like the conventional 2D disordered systems, there is a 
unique picture of the Anderson localization with four mo- 
bility edges in disordered graphene. We predict that An- 
derson MIT can be achieved in graphene by changing the 
carrier density to make a move of the Fermi surface across 
the mobility edges. In addition, we have found that the 
length scale of LDOS is considerably enhanced by the on- 
site interactions when the screening effects of interactions 
on the disorder potentials are also taken into considera- 
tion. 

Besides, the measure of Anderson localization by IPR 
has also been discussed, and a polynomial formula has 
been introduced for the lattice size scaling of IPR. As 
a result, precise localization length of the localized elec- 
tronic state can be obtained by the intercept found in the 
infinite-size limit. 
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Fig. 7. (color online) The disorder configuration and the 
corresponding LDOS are shown for the same parameters 
as in Fig. [6l except that the systems are all quarter- filled. 
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